The radial distributions of a heavy-light meson on a lattice 
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In an earlier work Q, the charge (vector) and matter (scalar) radial distributions of heavy- light mesons were 
measured in the quenched approximation on a 16^ x 24 lattice with (3 = 5.7, a lattice spacing of a ~ 0.17 fm, and 
a hopping parameter corresponding to a light quark mass about that of the strange quark. 

Several improvements are now made 

1. The configurations are generated using dynamical fermions with quark-gluon coupling /3 — 5.2 (a ~ 0.14 fm); 

2. Many more gauge configurations are included (78 compared with the earlier 20); 

3. The distributions at many off-axis, in addition to on-axis, points are measured; 

4. The data analysis is much more complete. In particular, distributions involving excited states are extracted. 

The exponential decay of the charge and matter distributions can be described by mesons of mass 0.9±0.1 and 
1.5±0.1 GeV respectively — values that are consistent with those of vector and scalar qg-states calculated directly 
with the same lattice parameters. 



1. INTRODUCTION 

A knowledge of the radial distributions of 
quarks inside hadrons could lead to a better 
understanding of the QCD description of these 
hadrons and possibly suggest forms for phe- 
nomenological models. As a step in this direction, 
we have measured the charge (vector) and matter 
(scalar) densities of a heavy-light meson. 

More explicitly, the heavy-light meson is sim- 
plified to being an infinitely heavy quark (Q) 
and an antiquark (q) with a mass approximately 
equal to that of the strange quark. The physi- 
cal meson nearest to this idealised meson is the 
Bs{5.37 GeV). We consider two flavours of dy- 
namical quarks {u and d). 

2. MEASUREMENTS 

2.1. Energies 

The basic quantity for evaluating the energies 
of a heavy-light meson is the 2-point correlation 



function - see Fig. |^ a) - 

C2(T) = ([/Q(x,t,r)P(x,t + r,t)), (1) 

where [/'^(x, i,T) is the heavy (infinite mass)- 
quark propagator and P(x, t + T,t) the light anti- 
quark propagator. The (...) means the average 
over the whole lattice. 

2.2. Heavy-light radial correlations 

When measuring the radial correlations, the 
basic quantity is the 3-point correlation function 
- see Fig. I b) - 

C3,e(™)(r,r) = (C/«PieP2), (2) 

where the Pi^2 sue the light anti-quark propaga- 
tors that go from the Q at time ti to the point r 
at i = and then return to Q at time — ^2- Here 
T = ti + 12- The probe at r is here considered 
to have two forms: i) G = 74 for measuring the 
charge distribution (c) of the q and ii) 8 = 1 for 
measuring the matter distribution (m). 
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Figure 1. Two- and three-point correlation func- 
tions C2(r) andC3,,(„)(T,r). 
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The strategy is to first fit the {C2{T)) by the 
approximate expression 



(72(T) = ^«"e- 



(3) 



Figure 2. The radial distribution of the ground 
state charge (C) and matter (M) densities. The 
lattice exponential fits are plotted with solid lines. 
Here r is in lattice units of a « 0.14 fm. 



This results in the energy eigenvalues and 
their eigenvectors (Ref. 

Given the and u", the C^^c{m){T, r) are then 
fitted by 

(53,cM(T,r) =^t;"e-^°*^a;"'5(r)e"^^*^i;^ (4) 

where the x"^{r) are varied giving directly the 
desired charge or matter density. In this work 
two levels of fuzzing with 2 and 8 iterations were 
used in addition to the unfuzzed lattice. This 
enabled us to extract information of the excited 
states as well. 

3. RESULTS 

3.1. Fits to the charge and matter densities 

It is convenient to parametrize the radial corre- 
lations in some simple way. We have considered 
pure exponential and Yukawa fits to the data as 
well as their lattice forms 



cos(r ■ q) 



(5) 



(lattice exponential) and 

TT ■r-^ cos(r • q) 



iL3 ^ 



(6) 



(lattice Yukawa). 

They all fit the data well, and the lattice forms 
are able to reproduce some of the structure (lat- 
tice artefacts) — see Fig. ^ The exponential de- 
cay of the charge and matter distributions can 
be described by mesons of mass 0.9 ±0.1 and 
1.5 ± 0.1 GeV respectively. The matrix elements 
with excited states are shown in Figs. ^, ^. 

3.2. Sum rules 

It is also of interest to consider the sum rules 
that sum over all values of r i.e. 



^c(m) 



(S.g3,c(^)(T,r)) 



(7) 



For the charge and matter densities this yields 
4 = 1.4(1) and /m = 0.9(1) respectively. The 
corresponding values for the quenched approxi- 
mation are Ic = 1.30(5) and /,n = 0.4(1). By 
charge conservation we should get /c = 1 in the 
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x^^/x^ ^ for the first excited state 



2 



-1 



\ c 






charge data 
» matter data 




M 


[ 









2 4 6 



r 



Figure 3. The ratio x^^/x^^ for charge (C) and 
matter (M) densities. The solid curves guide the 
eye. 



x^^/x" for the second excited state 
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continuum hmit. Renormahsation effects of f» 0.7 
enter due to the finite lattice spacing. 

4. CONCLUSIONS AND FUTURE 

The S'-wave charge and matter densities can 
be measured quite reliably out to » 0.8 fm and 
w 0.5 fm respectively. No difference can be 
seen when comparing dynamical and quenched 
fermions. This is probably due to the rather 
heavy mass of the light quarks. 

This work is by no means the last word on the 
subject, and a lot remains to be done. The next 
step could be, for example, measuring the P-, D- 
and F-wave densities. The correlations in the 
baryonic and {Q'^lf) systems should be studied as 
well. Wc would also like to understand the den- 
sities phenomenologically using the Dirac equa- 
tion. To check the continuum limit we should, of 
course, repeat the measurements using larger (3 
and bigger lattices. 
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Figure 4. The ratio x^^ /x^^ for charge (C) and 
matter (M) densities. 



